Multicomponent fluid of nonadditive hard spheres near a wall 
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A recently proposed rational-function approximation [Phys. Rev. E 84, 041201 (2011)] for the 
structural properties of nonadditive hard spheres is applied to evaluate analytically (in Laplace 
space) the local density profiles of multicomponent nonadditive hard-sphere mixtures near a planar 
nonadditive hard wall. The theory is assessed by comparison with NVT Monte Carlo simulations 
of binary mixtures with a size ratio 1:3 in three possible scenarios: a mixture with either positive 
or negative nonadditivity near an additive wall, an additive mixture with a nonadditive wall, and 
a nonadditive mixture with a nonadditive wall. It is observed that, while the theory tends to 
underestimate the local densities at contact (especially in the case of the big spheres) it captures 
very well the initial decay of the densities with increasing separation from the wall and the subsequent 
oscillations. 

PACS numbers: 61.20.Gy, 61.20.Ne,61.20. Ja, 68.08.De 



I. INTRODUCTION 

The study of mixtures near a fluid-solid interface is im- 
portant for the understanding of wetting and adsorption 
phenomena where competition among different compo- 
nents may occur. A simplified physical picture of ad- 
sorption may be obtained at a microscopic level if one 
considers the solid surface as a planar smooth hard wall 
confining the particles of the mixture. Thereby, one can 
describe the expected oscillations of the (partial) local 
particle densities in the neighborhood of the wall with an 
abundance of particles right at contact and a depletion 
nearby. Whereas confined fluid mixtures of additive hard 
spheres (AHS) have been widely studied within integral 
equation theories [H-Ill, Monte Carlo simulations [1-13|, 



and density- functional theories [1, IT3l - [23j . much less is 
known in the case of nonadditive hard spheres (NAHS) 

In a recent paper NAHS mixtures were stud- 

ied through the so-called rational-function approximation 
(RFA) technique [l^,!!^!, which amounts to choosing sim- 
ple (rational-function) expressions for the Laplace space 
representation of the radial distribution functions of the 
theory of liquids [Sll, • This allowed us to determine a 
nonperturbative, fully analytical (in Laplace space) ap- 
proximation. When the nonadditivity is set to zero, the 
approximation reduces to the Percus-Yevick (PY) ap- 
proximation for an AHS mixture. 

The purpose of the present work is to use the RFA 
scheme devised in Ref. [2^ to determine the structural 
properties of an n-component NAHS fluid near a hard 
wall interacting either additively or nonadditively with 
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the particles of the fluid mixture. A realization of the 
problem is obtained from a (n -f l)-component NAHS 
mixture, where one of the species, species 0, is taken to 
have a vanishing concentration and an infinite diameter. 
A similar approach was employed by Malijevsky et al. 
[isj to determine through the RFA the structural prop- 
erties of a multicomponent AHS fluid near an additive 
hard wall. In the present case, however, not only the 
particle-particle interaction may be nonadditive (i.e., the 
closest distance between the centers of two spheres of 
species i and j is in general different from the arithmetic 
mean of the respective diameters) , but also the particle- 
wall may be nonadditive as well. The latter possibility 
means that the closest distance from the planar wall to 
the center of a sphere may be different from the radius of 
the sphere. A similar problem has recently been consid- 
ered by Gonzalez et al. [s^ , where strong size selectivity 
is observed in a binary AHS mixture confined in a nar- 
row cylindrical pore such that each species of the mixture 
sees a different cylinder radius. 

We will compare our approximation results for the lo- 
cal density of particles at a distance z from the wall with 
exact canonical (fixed number of particles TV, volume V, 
and temperature T) Monte Carlo (MC) simulation re- 
sults for binary mixtures. In the simulation it is nec- 
essary to use two hard walls on the opposite far square 
faces of a parallelepiped simulation box with rectangu- 
lar lateral faces and to choose the two walls far enough 
so that bulk properties of the fluid can be extracted by 
looking at the center of the box. 

The agreement between theory and simulations is quite 
satisfactory. It is worse at contact (similarly to what hap- 
pens with the PY theory in the additive case but it 
rapidly improves as the distance from the wall increases, 
so that the flrst minimum (depletion region) and the sub- 
sequent oscillations are well predicted by our analytical 
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approach. To the best of our knowledge, our results con- 
stitute the first proposal for an analytical expression (in 
Laplace space) for the density profiles of a NAHS mix- 
ture confined by a (nonadditive or additive) hard wall. 
As such, the theory is expected to be useful to the exper- 
imentalist who needs easy formulas to determine profiles 
to compare with experimental data, thus bypassing the 
need of numerical experiments. 

The paper is organized as follows. In Sec. HI] we de- 
scribe the model of the confined fluid we are going to 
study. The RFA used to extract the structural prop- 
erties of the fluid is presented in Sec. IIIIl where some 
details of the wall limit are given in the Appendix. In 
Sec. IIVI we describe some details of the NVT MC simu- 
lation method we employed for confined binary mixtures. 
The results for the structural properties are presented in 
Sec. El where the RFA and our own MC simulation are 
compared. Finally, Scc. lVIl is left for concluding remarks. 



AHS (Ai2 = 0) or NAHS (A12 7^ 0) in the presence of a 
hard wall which, in addition, may interact additively or 
nonadditive ly with the fluid particles (see Sec. IIIIBj) . 



III. RATIONAL-FUNCTION APPROXIMATION 

A. General scheme 

In Ref. [13, the following proposal for the struc- 
tural properties of an n-component NAHS fluid dcflned 
through the Laplace transform Gij{s) of rgij{r) was 
given: 

n 

Gz,(s) =s-2^e-'^-^L,fc(s)Bfc,(s), (3.1) 



k=l 



with 



B-i(s) = l-A(s), 



(3.2) 



II. THE MODEL 

An n-component NAHS mixture in the d-dimensional 
Euclidean space is a fluid of Ni particles of species i with 
« = 1, 2, . . . , n, such that there are a total number of par- 
ticles Af = ^"^j^ A'i in a volume V, and the pair potential 
between a particle of species i and a particle of species j 
separated by a distance r is given by 



0.y(r) 



00, r < (Tij, 
0, r>(Ti,-, 



(2.1) 



where au = cr,; and (Ty = ^{ai + aj){l + A^), so that 
An = and Aij = Aji > —1. When Ay = for every 
pair i-j we recover the AHS system. In the present paper 
we will only consider the NAHS system in its single fluid 
phase. 

Let p = N/V be the total number density of the mix- 
ture and Xi = Ni/N be the mole fraction of species i. 
These are spatially averaged quantities that can differ 
from local values in conflned situations. 

The one-dimensional {d — 1) NAHS fluid admits an 
exact analytical solution for the structural and thermo- 
physical properties in the thermodynamic limit N ^ 00 
with p = N/V = const [M^. Moreover, the AHS fluid 
with d = odd allows for an analytical solution of the FY 
approximate theory [38l - l4l| . Such a solution in the case 
d = 1 reduces to the exact solution particularized to the 
additive mixture. 

Inspired by both the exact solution for one-dimensional 
NAHS mixtures and the FY solution for three- 
dimensional AHS mixtures, we have recently proposed 
an analytical approach for the three-dimensional NAHS 
system [2^. As said in Sec. HI the aim of the present paper 
is to use that approximation to determine the structural 
properties of a ternary mixture where one of the species 
(i ~ 0) is subject to the wall limit: xq — ^ and cfq — > 00. 
Such a ternary mixture represents a binary mixture of 



A,^.(s) = [A^,,(s)e"-^ - L,,(s)e-'^--] , (3.3) 



where I is the unit matrix. 



(3.4) 



L^s (1 - h,s) , 



1 



(di - CTj). 



(3.5) 
(3.6) 



Equations p.ip - (|3.5p provide the explicit s-depcndence 
of the Laplace transform Gy(s), but it still remains 

to obtain the two sets of parameters L^-^' and i^j^ 
This is done by enforcing the physical requirements 
\im.s^o s^Gij{s) = 1 and lims_i.os^^ [s^Gij(s) — l] = 
fH. The resuh is 



where 



1 - 7r/9*j 

(1 - npA,) (1 - TTp-^j) - TT^p^Pj\2.0^, 

npClj 

(1 - npAj) (1 - Trp^-j) - 7r2p2^j|2_o^ 



Aj = Pj\2,l ~ gMjIa.O; 



(3.7) 

(3.8) 

(3.9) 

(3.10) 
(3.11) 
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= ^^j\3,l ^ A^j|2,2 — ^^^^j\i,0^ (3-12) 



and we have called 



fc=l 



(3.13) 



As discussed in Ref. [28| , the inverse Laplace transform 
C^^ [Gij{s)] (r) may present a spurious behavior in the 
shell min((T.y, Ty) < r < max ( cr.y , Ty ) , where is the 
minimum of the list of values <Jik — a-kj {k = 1, . . . , n) that 



are different from aij. If (7^^ — a^j 



for all fc, then 



Tij — <Jij- The anomalous behavior of \Gij{s)] (r) for 
min(cry , Tij) < r < max((Ty , Ty) can be avoided with a 
series of corrections, the simplest one of which yields 



,-(r) = Q{r-a,j) 



C-'[G,A^)] (r) 



+a,e(7:y-r)(^-l)] , (3.14) 



where 



(3.15) 



Kij being the index associated with T27 , I.e ., Tjj — ^iKij 



The contact values are given by [2? 



.«(4) = f + o,(3i-i 



(3.16) 



The approximation p.l4p was referred to as RFA^^-* 



in Ref. 128] 

one has atk — akj = <Jij, so tnat Tij = 
gij{r) = r~^C~^ [Gij{s)]{r) coincides with the PY so- 
lution fHIIi. 



In the special case of AHS mixtures, 
dkj = Ci,-, so that Tii = an and 



B. Wall limit 

Now we assume that a single sphere of diameter ctq 
is introduced in the ?i-component fluid. This gives rise 
to an (n + l)-componcnt fluid, where the extra species 
(i — 0), being made of a single particle, has a vanishing 
concentration aio = in the thermodynamic limit 
00. With this proviso, Eq. (|3.ip can be easily extended 
to this {n + l)-componcnt mixture. 

According to Eq. p.3p . if xq = 0, the row i = of the 
matrix A is zero. As a consequence, the row i ~ and 
the column j = of the matrices and B have the 
forms 



Thus, application of Eq. p.ip to the pair i-0 with z > 1 
yields 

n 

G«(s) -s-^e— ««L,o(s) + ^G,j(5)^,o(s). (3.19) 

Therefore, the cross function Gio{s) (with i = 1, . . . , n), 
which is related to the spatial correlation between a par- 
ticle of species i > 1 and the single particle i = 0, is ex- 
pressed in terms of the matrix Gij{s) of the n-component 
mixture and the cross elements Lio(s) and Ajo{s). 

In principle, the nonadditivity of the i-0 interaction 
would be measured by the nonadditivity parameter A^o 
defined by = ^{o'o + o'i)(l + ^io)- However, the use 
of AiQ is not convenient in the wall limit (Tq — !■ 00 that 
we will take at the end. Instead, we define a nonad- 
ditivity distance Wi by ato = ^{o-q + (Ji) + Wi. Note 
that, since no 0-0 interaction is present, the definition 
of the diameter (Tq is somewhat arbitrary. In fact, if all 
Wi = w arc equal, the apparently nonadditive i-0 inter- 
action is indistinguishable from an additive interaction 
with CTo — > ctq -l- 2'w. Therefore, a true nonadditive i-0 
interaction requires, first, that n > 2 and, second, that 
not all {wi} are equal. Therefore, without loss of gener- 
ality, we take min(w,;; i = 1, . . . , n) = 0. This defines the 
diameter cfq unambiguously. 

As a next step toward the wall limit, we introduce the 
shifted radial distribution function 



= gioiz + cro/2). 



(3.20) 



Thus, while r is the distance between the centers of the 



pair i-0, 



■i(To represents the distance from the 



center of a sphere of species i to the surface of the single 
sphere j = 0. If we call ri(s) the Laplace transform of 
7i(z), the following relationship applies: 



G,o(s) 



yr,(.) 



(3.21) 



where r^(s) = dTi{s)/ds. 

Finally, we take the wall limit ag ^ 00. In that case, 
the function 7i(z) becomes the ratio between the local 
number density of particles of species j at a distance z 
from the wall, Pi{z), and the corresponding density in 
the bulk, p^. In an infinite system (as implicitly assumed 
in the theoretical approach) , the bulk and average values 
coincide, i.e., p° = pi(oo) = Xip. 

In the wall limit T[{s) can be neglected versus 
gqTi{s)/2 in Eq. dMD), so that 



r,(s)-2 lim a(7ie""^/2G,o(s) 



2Y,G,,{s)A,{s), 



Bia - -^lO, 



i > 1, 



B. 



iO 



fc=l 



BikAko, 



i > 1. 



(3.17) 



(3.18) 



(3.22) 

where in the second step we have made use of Eq. p.l9p 
and have defined 



Li{s) = lim Ctq L^is), 



(3.23) 



4 



Aj{s)= lim a^^e'^^'/^Ajois). (3.24) 

(To — >00 

These two quantities are evaluated in the Appendix. 

Once the Laplace transform ri{s) is well defined, let 
us consider the correction described by the second line 
of Eq. p.l4|) . First, we subtract ^(Tq to the distances, so 
that the shell min(tTio, Tjo) l£ r < max((Tio, Tio) becomes 
min (i(Ti + Wi,Ti) < z < max {^(Ji + Wi,Ti), where Ti is 
the minimum of the list of values aik — ^cr/c {k — 1, . . . ,n) 
that differ from jCTj + Wi. Again, Ti = \<Ji + Wi if <Jik — 
= ~\- Wi for all k. Finally, in the limit ~^ oo, 
one obtains 

7.(^) = e (^2 - ia, - [C-^ {T,{s)] (r) 

+2ae(f, -z)(n-z)}, (3.25) 

with 

(3.26) 

where Ki is the index associated with Ti, i.e., r,; = aiK^ — 
ifTKi, s-iid the quantities L'^"^ and L.^^-* are defined in the 
Appendix. 

The inverse Laplace transform in Eq. p.25p can be 
easily performed numerically p^ . On the other hand, 
the density ratio ji at the shortest distance from the 
wall z — ^ai + Wi can be derived analytically. From Eq. 
(|3.16p we easily obtain 

7i = ^o-^ + w'i^ = 2Zf ^ + 2Ci - ifTi - . 

(3.27) 

The fact that the general scheme gives well defined ex- 
pressions in the wall limit {xq = 0, ctq — > oo) is a stringent 
test on the internal consistency of the RFA approach. It 
also shows the convenience of dealing with explicit, ana- 
lytical expressions from which the subsequent limits can 
be taken. 



IV. MONTE CARLO SIMULATIONS 

We have simulated a binary mixture {n ~ 2) of NAHS 
through canonical NVT MC simulations in a box of fixed 
volume and sides L^, Ly, and Lz with = Ly and 
Lz S> Lx- Periodic boundary conditions are enforced 
along the x and y directions, but two impenetrable hard 
walls are located at z = and z ~ L^. The particles 
are initially placed on a simple cubic regular configura- 
tion along the z direction with a first crystal layer of 
particles of species 1 juxtaposed to a crystal layer of 
particles of species 2. We reject the ith particle move 
only in case of overlap with any other particle, i.e., if 
Tij < o'ij for some j, or with one of the walls, i.e., if 
niin{zi, Lz — Zi) < ^ai + Wi. The system is then equili- 
brated for 10^ MC steps (where a MC step corresponds 



to a single particle move) and the properties are gen- 
erally averaged over additional 10^ MC steps for pro- 
duction. The maximum particle displacement, the same 
along each direction, is determined during the first stage 
of the equilibration run in such a way as to ensure an 
average acceptance ratio of 50% at production time. As 
a compromise between the condition Lz 3> L^ and the 
computational need of not having too high a number 
of particles, we have taken L^ ~ 10<ti and Lz = 30cri 
in all the simulations presented, except a control case 
with Lz = 60cti (see below). The local density profiles 
7i(z) = pi{z)/ p° are obtained, for each species, from his- 
tograms of the z coordinates of the particles in bins of 
width 0.01 (71. The bulk values are evaluated in the 
region of the simulation box with z « Lz/2, where a neg- 
ligible infiuence from the walls is expected. Due to the 
finite value of Lz, the bulk total density p^ ^ p\ + p^ 
and the bulk mole fraction x\ = P\l P° differ from their 
respective average values p and x\. 

V. RESULTS 
A. Representative systems 

In the binary case, there are five independent dimen- 
sionless parameters of the problem: the size ratio o^j <J\, 
the particle-particle nonadditivity parameter A12, the 
particle- wall nonadditivity parameter max(u;i, ii;2)/CTi 
[remember that, by convention, T[im{w\,W'i) = 0], the 
average mole fraction x\, and the average reduced den- 
sity pa\. Here, a\ is chosen as the diameter of the small 
spheres and henceforth it will be used to define the length 
unit. 

In order to focus on the nonadditivity parameters, we 
have chosen Uijax = 3 for all the systems. Next, three 
classes of systems have been considered: (i) a nonaddi- 
tive mixture in the presence of an additive wall, (ii) an 
additive mixture in the presence of a nonadditive wall, 
and (iii) a nonadditive mixture with a nonadditive wall. 
As representative examples of class (i) we have chosen an 
equimolar mixture with either positive (system A of Ta- 
ble |I| or negative (system B of Table |I| nonadditivity and 
a mixture with an excess of small spheres and negative 
nonadditivity at two densities (systems CI and C2, re- 
spectively). As examples of class (ii), we have chosen an 
equimolar mixture where the wall presents an extra re- 
pulsion to either the large spheres (system D) or the small 
spheres (system E). Finally, class (iii) is represented by 
system F, which is analogous to system D, except that 
the mixture has a negative nonadditivity. The reduced 
densities pcrf range from to ^, so that the total num- 
ber of particles N = pL^Lz ranges from 100 to 600. It 
is also convenient to measure the density in terms of the 
effective packing fraction r^eff = ^pJ2i 

related 

to van der Waals's one-fluid theory [i^, whose values are 
indicated in the last column of Table ID In the low-density 
regime, two mixtures with the same value of ^off would 
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TABLE I. Values of the nonadditivity parameters (A12, wi, 
and W2), the average mole fraction (2:1), and the average den- 
sity (p) for the representative systems considered in this work. 
In all the cases 02 1 ox — 3. The table also includes the values 
{x\ and p^) measured in the bulk region z ~ Lz/2 in our MC 
simulations with Lz/o"i = 30. 



Label 


A12 


wi/ai 


102/ ai 


Xl 


4 


pal 




^cff 


A 


0.2 








0.5 


0.469 


1 

30 


0.0337 


0.243 


B 


-0.2 








0.5 


0.503 


1 

20 


0.0513 


0.237 


CI 


-0.2 








0.9 


0.896 


1 

10 


0.1025 


0.095 


C2 


-0.2 








0.9 


0.898 


I 
5 


0.2040 


0.190 


D 








0.35 


0.5 


0.475 


1 

30 


0.0345 


0.192 


E 





0.35 





0.5 


0.511 


1 

20 


0.0503 


0.288 


F 


-0.2 





0.35 


0.5 


0.486 


1 

30 


0.0350 


0.158 
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012345678 



have the same compressibility factor. 

B. Bulk values 

The bulk values x\ and measured in the MC sim- 
ulations with = 30(71 are also included in Table HI In 
all the cases the bulk density is larger than the aver- 
age density p. This is due to the fact that the effective 
length available to the spheres of species i is not Lz but 
Lz — (ci + 2wi). As a consequence, the larger deviation 
between p^ and p takes place for systems D (3.5%) and 
F (5.0%), i.e., those systems where the walls produce an 
extra repulsion {w2/(7i = 0.35) on the big spheres. This 
compression effect is only partially compensated by the 
accumulation of particles at contact with the walls. In 
the case of the bulk mole fraction x\ , the situation is less 
obvious. Note the identity 

4S-7^r^' ^-7i(^). (5.1) 

Even if the right-hand side of Eq. ()5.ip is generally 
smaller than 1, the fact that p^ > p can give rise to 
x\ < xi; i.e., the bulk would be richer in big spheres than 
on average. This is what actually happens for systems 
A, CI, C2, D, and F. This effect is especially important 
in systems A and D since in those cases the right-hand 
side of Eq. (|5.ip turns out to be larger than 1 (see Figs. 
[1] and [5] below for a visual confirmation). Exceptions to 
the property Xi < xi arc represented by systems B and 
E. In those cases, the right-hand side of Eq. (|5.ip is suf- 
ficiently smaller than 1 (see Figs. [3] and IH] below) as to 
compensate for the ratio p° lp> 1. 

Now we turn our attention to the density profiles. 
When presenting the theoretical RFA results for each sys- 
tem we have used two criteria. In the first criterion, the 
quantities p and x\ appearing in the theoretical scheme 
described in Sec. lIIII have been identified with the average 
values employed in the simulations. In the second crite- 
rion, the RFA quantities p and x\ have been identified 



FIG. 1. (Color online) Local density profiles 7^(2) = pi{z)/ fhl 
for system A {a2la\ = 3, A12 = 0.2, wi = UI2 = 0, X\ = 
0.5, pa\ — 1/30). The lines represent the RFA theoretical 
predictions using the average values x\ and p (dashed lines) or 
the empirical bulk values x^ and p° (solid lines). The symbols 
represent our MC simulations with L^jax = 30 (circles) or 
Lz/ai = 60 (squares). The inset shows the MC data in the 
whole domain Q < z < with Lz/a\ — 30. In the MC 
results, the error bars are within the size of the symbols used 
in the graph. 

with the bulk values p^ and x\ found in the MC simu- 
lations with Lz = 30(Ti. As said before, the theoretical 
approach deals with formally infinite systems {Lz 00) 
where the average and bulk quantities coincide. However, 
when making contact with simulation data correspond- 
ing to finite Lz the use of cither the average or the bulk 
values in the RFA may be important. 



C. Nonadditive mixture and additive wall 

Figure [1] shows the MC and RFA results for the two 
(relative) density profiles 7^(2) = pi{z)/p° {i = 1,2) in 
the case of system A (positive nonadditivity). In this 
system 77 = 0.9tTi > ifTi and T2 = 1.9(Ti > ^(T2, so that 
the correction term given by the second line of Eq. p.25p 
is used in the RFA curves. 

The inset of Fig. [1] shows the MC results for both den- 
sity profiles in the whole region Q < z < Lz ~ 30(Ti. We 
can see that the separation between both hard walls is 
large enough as to identify a well defined bulk region in 
the center. We have chosen system A to assess the in- 
fluence of flnite Lz by carrying out a control simulation 
with Lz = 60(71. The new bulk values are x\ = 0.485 and 
p'Vi = 0.0334, which, as expected, are closer to the av- 
erage values than in the case Lz = SOcti (see Table |l|. As 
seen from Fig. [U the MC data obtained with Lz = SOcti 
and Lz = 60cti are hardly distinguishable, except near 
contact where the smaller system, having a larger bulk 
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FIG. 2. DifTerences A7i(z) = 7i(z)l3o — 7«(-^)l60 between the 
local densities "/iiz) = pi{z)/p° evaluated in MC simulations 
with L^/ai = 30 and those evaluated in MC simulations with 
Lz/(Ji = 60 for system A (02! o\ = 3, A12 = 0.2, w\ = W2 — 
0, si = 0.5, pal = 1/30). 



density, presents slightly higher values of ^fi(z). A more 
detailed comparison is made in Fig. [51 where the dif- 
ferences between the values of 7^(2:) as obtained with 
both values of are shown. Figure [2] confirms that the 
smaller system (L^ = 30(Ti) presents larger values for the 
two reduced densities near contact than the larger sys- 
tem (L^ = 60(Ti). For higher separations the differences 
are much less important, but yet it is interesting to note 
that the smaller system tends to present larger values of 
72(2:) but smaller values of 71(2;). 

Now let us go back to Fig. [T] and comment on the per- 
formance of the RFA. We observe that the RFA under- 
estimates the local densities at contact (i.e., at z = \(Ji)- 
On the other hand, the decay of the local densities near 
the walls and the subsequent oscillations are very well 
captured by the theory. It is interesting to remark that 
the agreement with the MC data near contact improves 
when the bulk values instead of the average ones are used 
in the theory. 

The profiles for system B (negative nonadditivity) 
are displayed in Fig. |3l In this case 17 = O.Icti and 
T2 = l.lci. Since < \<Ji, the correction term in the 
second line of Eq. (|3.25p vanishes. Comparison between 
Figs. [Hand [3] shows that, in going from system A to sys- 
tem B, the local variation of the density of the big spheres 
is enhanced, while the local density of the small spheres 
becomes less structured. Here there are two competing 
effects at play. On the one hand, at a fixed density, the 
change from positive to negative nonadditivity produces 
a weaker density structure near the wall, as the exact re- 
sult to first order in density clearly shows. On the other 
hand, at a fixed nonadditivity, an increase in density in- 
duces a higher structure. It seems that, in the transition 
from system A to system B, the latter effect dominates 



FIG. 3. (Color online) Local density profiles 7^(2) = pi{z)/ fhl 
for system B [02/01 = 3, A12 = —0.2, tui = W2 = 0, Xi = 
0.5, peri = 1/20). The lines represent the RFA theoretical 
predictions using the average values xi and p (dashed lines) or 
the empirical bulk values x\ and p° (solid lines). The symbols 
represent our MC simulations with L~/g\ = 30. In the MC 
results, the error bars are within the size of the symbols used 
in the graph. 



in the case of the big spheres (which arc very weakly 
influenced by the small component) and the former ef- 
fect docs it in the case of the small spheres (which are 
strongly influenced by the presence of the large compo- 
nent). It is interesting to note that all these features are 
very well described by the RFA, especially in the case of 
71(2:). The contact value of 72 is better estimated in sys- 
tem A than in system B, while the opposite happens for 
the contact value of 71. Note also that a small discrep- 
ancy is observed near the second peak of 72(2) in Fig. [31 
For this system the RFA is practically insensitive to the 
use of the bulk values instead of the average ones. 

In systems A and B the big spheres occupy as much 
as 27 times more volume than the small ones, so the 
global properties of the mixture are dominated by species 
2. A more balanced situation takes place in systems CI 
and C2, where the ratio of partial packing fractions is 
X2(J2/ xiai ~ 3. In these cases the high concentration 
asymmetry requires a long simulation run time to reach 
thermal equilibrium for the big spheres. 

The results for systems CI and C2 are shown in Fig. 
[H At the smaller density (system CI) the agreement 
between theory and simulation is almost perfect. As the 
density is doubled (system C2), some small deviations are 
visible, especially in the case of the big spheres. Again, 
the RFA with the bulk values behaves near contact better 
than with the average values. 
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FIG. 4. (Color online) Local density profiles 71(2) = pi{z)/ p° 
for (a) system CI {02/01 — 3, A12 = —0.2, wi — 'W2 = 0, 
xi_ = 0.9, pal = 1/10) and (b) system C2 (0-2/cri = 3, A12 = 
—0.2, ^1=^2 = 0, xi = 0.9, pcrf = 1/5). The lines represent 
the RFA theoretical predictions using the average values x\ 
and p (dashed lines) or the empirical bulk values x\ and p^ 
(solid lines). The symbols represent our MC simulations with 
Lz/ai = 30. In the MC results, the error bars are within the 
size of the symbols used in the graph. 



FIG. 5. (Color online) Local density profiles 7^(2) = pi{z)/ fhl 
for system D (a^jox = 3, A12 = 0, wi = 0, 1112/01 = 0.35, 
xi — 0.5, pa'l = 1/30). The lines represent the RFA theo- 
retical predictions using the average values x\ and p (dashed 
lines) or the empirical bulk values x\ and p° (solid lines). 
The symbols represent our MC simulations with L^/ai = 30. 
In the MC results, the error bars are within the size of the 
symbols used in the graph. 



thus, the contact value and the oscillations are more pro- 
nounced in system E. These effects are enhanced by the 
larger density of system E relative to that of system D. 
However, 71(2:) near contact is higher in system D than 
in system E, so that the effect of wall nonadditivity com- 
pensates for the increase of density in the case of the 
small spheres, analogously to what happens with sys- 
tems A and B (see Figs. [T] and [3]). All these features are 
correctly accounted for by the RFA, although the quan- 
titative agreement near contact is again worse than that 
after the first minimum, especially in the case of 72 (2;) ■ 
Note also that the influence on the RFA curves of the 
use of the bulk versus the average values is noticeable in 
system D but not in system E. 



D. Additive mixture and nonadditive wall 

Now we consider the cases where the mixture is ad- 
ditive but the wall treats differently both species. The 
extra repulsion affects the big spheres in system D and 
the small spheres in system E. In both cases Ti < ^ai+Wi, 
so that again the correction term in Eq. p.25|) does not 
apply. 

The results for systems D and E are shown in Figs. 
[5] and [SI respectively. In the case of system D there is 
much more room for the small spheres to sit between the 
wall and the big spheres than in the case of system E. 
As a consequence, the big spheres "feel" the presence of 
the wall more in the latter case than in the former and. 



E. Nonadditive mixture and nonadditive wall 

The more general situations where both the particle- 
particle and the wall-particle interactions are nonadditive 
is, of course, richer than the preceding classes. As a sim- 
ple representative system we consider the same case as 
in system D (wall additionally repelling the big spheres), 
except that, in addition, species 1 and 2 interact with 
negative nonadditivity. The resulting system F (see Ta- 
ble |l| is also close to system B, except that now the wall 
is nonadditive and the density is smaller. As in systems 
B-E, the correction term in Eq. p.25p is not needed. 

The local densities for system F are plotted in Fig. 
[T] Comparison with Fig. [S] shows that the density pro- 
file of the big spheres is practically unaffected by the 
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FIG. 6. (Color online) Local density profiles 7; (2) = pi(z)/ 
for system E {02/01 — 3, A12 — 0, wi/oi = 0.35, W2 ~ 0, 
= 0.5, po'i = 1/20). The lines represent the RFA theo- 
retical predictions using the average values x\ and p (dashed 
lines) or the empirical bulk values x\ and p° (solid lines). 
The symbols represent our MC simulations with L~/o\ — 30. 
In the MC results, the error bars are within the size of the 
symbols used in the graph. 



FIG. 7. (Color online) Local density profiles jiiz) = pi{z)/ fhl 
for system F {02 /oi = 3, A12 = —0.2, wi = 0, W2/01 = 0.35, 
x\ — 0.5, poi = 1/30). The lines represent the RFA theo- 
retical predictions using the average values x\ and p (dashed 
lines) or the empirical bulk values x\ and p° (solid lines). 
The symbols represent our MC simulations with L^/oi = 30. 
In the MC results, the error bars are within the size of the 
symbols used in the graph. 



nonadditive character of the 1-2 interaction. This is not 
surprising taking into account that, as said before, the 
big spheres occupy 27 times more volume than the small 
ones and, therefore, the presence of the latter has little 
impact on the properties of the former. On the contrary, 
the nonadditivity has a large influence on the local den- 
sity profile 71 (z). Since spheres of species 1 and 2 can 
overlap to a certain degree in system F, the big spheres 
partially alleviate the influence of the wall on the small 
spheres with respect to the case of system D. As a con- 
sequence, the local density of the small spheres is less 
structured in system F than in system D. Like in system 
D, the RFA performs very well in system F, especially 
when the bulk values are used. 



VI. CONCLUSIONS 

In this work we have developed a simple analytical (in 
Laplace space) nonperturbative theory for the local den- 
sity profiles of a multicomponent fluid of NAHS confined 
by an additive or nonadditive hard wall. The theoretical 
approach is based on the specialization of the RFA tech- 
nique recently proposed [28| to the case where an extra 
single particle of diameter cto is added to the mixture and 
then the limit of an infinite diameter uq ^ 00 \s taken. 
The RFA reduces to the exact solution of the FY ap- 
proximation for zero nonadditivity, both in the particle- 
particle and in the particle- wall interactions, but remains 
analytical even when nonadditivity prevents one from ob- 
taining an analytical solution of the PY theory. 



While the theory applies to any number of components, 
we have focused on a binary mixture with a size ratio 1 : 3 
plus a hard wall. This has allowed us to compare the 
theoretical results against exact MC simulation. Several 
representative scenarios have been considered (see Table 
a positive (system A) or negative (systems B, CI, and 
C2) NAHS fiuid with an additive wall, an AHS mixture 
with a nonadditive wall pushing either the big (system 
D) or the small (system E) spheres, and a NAHS mix- 
ture with a nonadditive wall (system F). In all the cases, 
a reasonably good agreement between our theory and the 
MC simulations have been found for the (relative) par- 
tial local densities ^i{z). The agreement is worse near 
contact, where the RFA underestimates the MC values, 
but rapidly tends to improve for larger distances, so that 
the initial decay of the local densities and the subsequent 
oscillations are rather well captured. Note that, since the 
RFA can be seen as a sort of continuation of the AHS PY 
solution to the NAHS realm [2^ , it is not surprising that 
some of the features of the PY solution remain. One of 
those features is the underestimation of the contact val- 
ues Another PY feature, namely the possibility of 
predicting a negative first minimum at sufficiently high 
densities, is also inherited by the RFA. 

As shown by Figs. [T]and[3H71 the performance of the 
RFA is usually better in the case of the small spheres 
{i = \) than for the large spheres {i = 2). This is in 
part due to the physical observation that the local den- 
sity structure of species 1 is milder than that of species 
2. Another technical reason has to do with the fact that, 
while the separation between both walls is sufficiently 
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large for the spheres near a wall not to be much influenced 
by the presence of the other wall, the unavoidable "com- 
pression" effect is more important for the big spheres 
{L2/U2 = 10) than for the small spheres {Lz/cfi = 30). 
As Fig. [2] illustrates, when the separation between both 
walls is doubled, the effect on the density near the walls 
is more pronounced for the big spheres than for the small 
ones. Finite-size effects are also related to the small dif- 
ferences between the average densities and their hulk val- 
ues in the central region z ~ L^/'i. We have checked that 
our theoretical approach exhibits a slightly better agree- 
ment with simulations when the empirical bulk values are 
used instead of the average values. 

Our theory, being a simple analytical one, can be ef- 
ficiently used to easily extract many-body approximate 
properties for confined fiuids under other interesting situ- 
ations different from the representative ones examined in 
this work. For instance, extreme cases like the Widom- 
Rowlinson jSjicj {a i = (T2 = with 012 finite) or the 
Asakura-Oosawa [13,1111 (ci = and A12 > 0) confined 
fiuids can be studied. Another avenue for application 
of the RFA is the depletion potential between two big 
spheres immersed in a sea of small spheres (49} interact- 
ing nonadditively with them. 



so that the denominator becomes 



D=[l- -pAh.o 



p2 (Ml^ - M2,oA/2,2) . (A.5) 



Equations (IA.ll) - (jA.5|) apply to any value of ctq. From 
Eqs. (|3J | - ((3^ it is easy to sec that both ij^' and ijg' 
are linear functions of ctq. Thus, taking the limit p.23p 
one gets 



with 



2,0 



2D 



(A.6) 
(A.7) 



L. 



(1) 



2D 



1 - (M3,o + 3Af2,i) + (y + T^pM2,0 

(A., 



Analogously, from Eqs. ([53]) and ([53)1 . the limit ([5^ 
becomes 



Appendix: Evaluation of Li{s) and Aj{s) 



Us) = 



iVj(s) -ij(s)e-('"^+'"^)"l , (A.9) 



Let us recall that = \{(Ti + ctq) + Wi with wi > 0. 
Therefore, according to Eq. (|3.6p . feio — + Wi. Thus, 
Eqs. (i3J0l) - (j3T2| yield 

Ao = yM2,o + ^M3,o + ^M2,i, (A.l) 



where 



1 - {aj +Wj)s + 



(A.IO) 



'i'0 = -yM2,0 + ^Af3,0-^M2,i, (A.2) 
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k=l 

Interestingly enough, the terms proportional to co a-nd 
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